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ABSTRACT 

We present an observational study of the protostehar core B335 harboring a low- 
mass Class source. The observations of the H^^CO^(J = 1-0) line emission were 
carried out using the Nobeyama 45 m telescope and Nobeyama Millimeter Array. Our 
combined image of the interferometer and single-dish data depicts detailed structures of 
the dense envelope within the core. We found that the core has a radial density profile 
of n(r) oc r"'' and a reliable difference in the power-law indices between the outer and 
inner regions of the core: p ~ 2 for r > 4000 AU and p ~ 1.5 for r < 4000 AU. The dense 
core shows a slight overall velocity gradient of ~ 1.0 km s^^ over the scale of 20,000 AU 
across the outfiow axis. We believe that this velocity gradient represents a solid-body- 
like rotation of the core. The dense envelope has a quite symmetrical velocity structure 
with a remarkable line broadening toward the core center, which is especially prominent 
in the position-velocity diagram across the outflow axis. The model calculations of 
position-velocity diagrams do a good job of reproducing observational results using the 
collapse model of an isothermal sphere in which the core has an inner free-fall region 
and an outer region conserving the conditions at the formation stage of a central stellar 
object. We derived a central stellar mass of ~ O.IMq, and suggest a small inward 
velocity, ■Wr^nnf ~ Okms~^ in the outer core at > 4000 AU. We concluded that our data 
can be well explained by gravitational collapse with a quasi-static initial condition, such 
as Shu's model, or by the isothermal collapse of a marginally critical Bonnor-Ebert 
sphere. 

Subject headings: ISM: clouds — ISM: individual objects (B335) — ISM: molecules — 
stars: formation 
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1. INTRODUCTION 

In order to understand the formation processes of low-mass stars, it is important to investigate 
the properties of dense (~ lO^cm"'^) cores in molecular clouds. Such compact (~ 0.1 pc) cores 
supply material to newly forming stars through dynamical gravitational collapse, however, detailed 
physical processes are still uncertain. One of the investigative approaches is to derive the detailed 
density and velocity structures from observations of (pre-)protostellar cores which are expected to 
retain more inforrnation than Class I /II objects for the initial conditions of gravitational collapse 
(|Andre et al.lll99.4 ISaito et al.l[l99d : IPuruva et allboOfil ^. 



Dust continuum emission imaging at millimeter and submillimeter wavelengths using single- 
dish radio telescopes has revealed the radial density profiles, p(r), of (pre-)protostellar cores 



(jWard-Thompson et al.lll994 



199S 



Andre et al 



1996 



Shirlev et al.l 120001 ). Recent investigations 



have demonstrated that the pro files of Class O/I so ur ces can be fit t ed by single power-law profiles 
over a wide range of radii (e.g., IShirley et al.ll2000l ). IShirley et al.l ([2002) modeled Class source 
maps using a single power-law density distribution p{r) oc and found that most of them can 
be well fitted with a power-law index of p ~ 1.8. Cas kinematics in dense cores have been investi- 
gated by molecular line observations. It has been shown that asymmetric double-peaked profiles of 
optically thick lines detected toward star-forming cores are considered to be a signature of collapse 



Lee et al. 


(1999. 


2001 


2004) 



simple two-layer radiative transfer model and suggested an inward motion of subsonic speed (0.02- 
O.lkms"^) that extended to 0.1 pc in the pre-protostellar core L1544. The velocity structures of 
dense cores, including not only infalling motion but also rotation, have also been investigated using 
first and se cond moment ima ges of line emission taken with single-dish telescopes and interferom- 
eters (e.g., 



Tobin et al 



2011 



Statistical studies for the rotation of dense cores were conducted 



using NII3 and N2II"'" line emissions by iGoodman et al.l (119931 ) and ICaselli et al.l (120021). in which 



the ty pical velocity gradients are found to be l-2kms~^ pc~^. On the other hand, iChen et al 
(j2007l ) observed N2H"*"(J = 1-0) line emission toward nine low-mass protostellar envelopes down 
to ~ lOOOAU scales. The mean velocity gradient estimated in t heir samples is_~ 7kms~^ pc^^, 
which is much larger than the velocity gradients of dense cores. iTobin et al.l (1201 ll ) analyzed the 
kinematics of 17 protostellar systems. They found that the velocity gradients obtained with inter- 
ferometric data (~ 8.6kms~^ pc^^) are considerably larger than those that also have single-dish 
data (~ 2.3 kms~^ pc~^), which indicates accelerating infall and spinning-up rotational velocities 
toward the core center. 

The environments surrounding forming stars are composed of structures with different scales: 
circumstellar disks (< 100 AU), infalling envelopes (~ 1000 AU), and dense cores (~ 10,000 AU). 
Recently, an approach that combines data obtained with single-dish telescopes and interferometers 



has been widely used to investigate th e physics i n pro tostellar systems (e.g., iFuruva et al 



2006 



Takakuwa et al.ll2007l : lYen et al.ll201lh . lYen et al.l (|201lh performed C^«0(J = 2-1) and CS(J = 7 



6) observations toward the protostellar envelope of B335 with the Submillimeter Array and single- 
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dish telescopes, and imaged by combining those data. They derived the specific angular momentum 
of the envelope and found that specific angular momenta tend to be larger as evolution progresses 
by comparing with other Class 0, I, and II sources. 

Theoretically, two extreme models for the core evolution have been proposed for isolated low- 
mass star formation. T he similarity solution of Lar son-Penston describes the density evolution of 
isothermal gas spheres (jLarsonlll969l : |Penstonlll969l ). When a central object is formed (t = 0), the 
gas sphere reaches the density profile of p{r) = 4.4(Cg/27rG)r~^ and velocity field of v{r) = 3.3 Cg, 
where Cg is the isothermal sound speed and G is the gravitational constant. This model of the 
Larson-Penston solution is refer red to as " runaway" collapse. On the other hand, the isothermal 
similarity solution proposed by IShul (| 19771 ) describes a core that is slowly increasing its central 
density through ambipolar diffusion while maintaining kinematic balance (i.e., v{r) = 0), and 
moving toward dynamical collapse. The density profile achieves a singular isothermal sphere, 
p{r) = (Cg/27rG)r~^, at t = 0, which is the initial condition of dynamical collapse after protostar 
formation. This model is the most static and is referred to as "inside-out" collapse. M oreover, 
extens ions for t > of the Larson-Penstq n solution and generalization were developed by iHunter 
(jl977l ) and IWhitworth &: SummersI (jl985l ). At later times (t > 0), in both of the solutions, the 
density and velocity structures attain a free-fall profile, p oc r~^/^ and v oc r~^/^, respectively, from 
the center to the outside, with a sound speed for the Shu solution that accompanies the rarefaction 
wave and with a supersonic velocity for the Larson-Penston solution. The mass infall rate is 
predicted to be ~ c^/G for the Shu solution and 48 times higher than this for the Larson-Penston 
solution. 

In this paper, we present an observational study of the dense core associated with a Class 
protostar within B335. The Bok Globu le B335, otherwise known in the literature as CB199 in the 
catalog of IClemens &: Barvainia (|l988l ) or L663, is a typical low-mass star-forming region. B335 
appears as an opaque dark cloud on optical images and is one of the best candidates for study- 
ing the initial conditions of star formation, because it is isolated fr om other star-for ming regions 
and is near the Sun. In this paper, we adopt a distance of 150 pc (IStutz et al.ll2008l ). B335 con- 
tains a far-infrared (F IR) source, IRAS 19347-1-0727, which is bright in submillimeter wavelengths 
(jChandler et al.lll99d ) and shows a combination of characteristics that indicates one of the clearest 
exam ples of very young st ars. This FIR source is associated \yith the dense molecula r gas enve- 



2010 



lope ( Frerking et al.l 119871 : iMenten et al.l Il989l : iHasegawa et al.l Il99ll : ISaito et al.l Il999l : lYen et al 



201 ih. and is cl assified as a Class protostellar object based on its spectral energy distri- 



bution (lBarsonvlll994l'). This source is considered to be the driving source of the bipolar ou tflow 
(jPrerking Langeilll98 j : Iffirano et al]ll988l . [l993 : Ichandler k Sargentlll993l : IYcu et ahlboiol l. ex- 
tending from east to west with an inclination angle from the sky plane of ^ 10° and an opening 
angle of ~ 45°. Toward the IRAS source in B335, IZhou et al.l (jl993l . Il994l ) have obtained blue- 
skewed line profiles in the CS(J = 2-1, 3-2, and 5-4) and H2C0(J = 2i2-lii and 3i2i:2ii_)_Jine 
emissions, which can be exp lained by a spherically symmetric inside-out collapse model ( Shulll977l ). 



Moreover, IChoi et al.l (jl995l ) computed radiative transfer using a Monte Carlo method with a model 
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of inside-out collapse. They confirmed that the double-peaked line profiles indicate an infa ll mo- 
tion and derived the physical parameters of the collapsing envelope. IVelusamv et al.l (jl995l ) made 
interferometric observations in the CCS line and also interpr eted the kinematics as an infalling 
motion, whereas the interferometric images of CS(J = 5-4) bv lWilner et al.l (|200d ) were found to 
be dominated by clumps associated with the outflow cavity. 

As mentioned above, the detailed kinematics of infall motion provide crucial information about 
the properties of the gravitational collapse of the core, which we can compare with the theoretical 
pictures and distinguish between them. Few observational studies until now, however, have been 
capable of unveiling the velocity structures within dense cores in the early phases of star formation 
over a wide range of spatial scales. We carried out the observations with the Nobeyama 45 m 
telescope and the Nobeyama Millimeter Array (NMA) in the H^'^CO^(J = 1-0) line emission. 
Observations with the 45 m telescope trace the overall density and kinematic structures of the 
star-forming core. For the high-resolution observations with NMA, we combined the data with the 
45 m telescope data to image the structures of the dense core over a wider spatial frequency range 
down to the inner dense envelope scales. 



OBSERVATIONS AND DATA ANALYSIS 



2.1. Nobeyama 45 m Telescope 



The mapping observations of the H^^CO+(J = 1-0) line emission at 86.75433 GHz toward the 
B335 regi on were carried out during 2006 March-May . We used the 25 Beam Array Receiver System 
(BEARS: ISunada et al.ll2000l : lYamaguchi et al.ll200d ). which consists of superconductor-insulator- 
superconductor (SIS) receivers on a 5 x 5 grid with a separation of 41'.'1, in the double sideband 
(DSB) mode. The observations were made in the Qn-The -Fly (OTF) mode of the Nobeyama Radio 
Observatory (NRO) 45 m telescope (ISawada et al.l 120081 ). The mapping center was placed at the 
IRAS source 19347-^0727 associated with B335, R.A. = 19'^34°>35!1, deck = -h07°27'20'.'0 (B1950), 
and the mapping region was ~ 4 x 4arcmin^. We obtained five OTF-scan data over this region in 
the R.A. and deck directions in each and merged them to create a final map cube. At 87 GHz, the 
half-power beam width and main-beam efficiency were 18'.'5 and 0.5, respectively. As the back end, 
we used 25 sets of 1024 channel autocorrelators (AGs), with a frequency resolution of 31.25kHz 
corresponding to a velocity resolution of 0.108kms~^. The system noise temperatures were between 
200 K and 280 K. We took the off position of R.A. = 19''36™35!2, deck = -F07°27'22'.'0 (B1950) and 
used the standard chopper wheel method to convert the receiver output into intensity scale. 
The telescope pointing was checked once an hour by five-point observations of the SiO maser from 
RT Aql in the 43 GHz using the SIS receiver (S40). The pointing accuracies were within 3" for the 
H^^GO"*" observations. 
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For the H^^CO"'' data obtained with the OTF observations, IDlJj-based reduction software, 
NOSTAR (Nobeyama OTF Software Tools for Analysis and Reduction), was used for flagging, 
baseline subtraction, and making the map cube. We corrected for relative gain differences among 
the 25 beams using correction factors provided by NRO, and then made the baseline fitting and 
subtraction. We made a three-dimensional image cube that has an effective spatial resolution of 
25'.' 1 and an achieved noise level in channel images (la) of 95 mK in with a velocity resolution 
of 0.108 kms-^ 

We also observed three inversion spectra of NH3(J, JC) = (1,1), (2,2), and (3,3) at 23 GHz 
toward the IRAS source associated with B335 using the HEMT receiver (H22) during 2006 March- 
May. At 23 GHz, the half-power beam width and main-beam efficiency were 78" and 0.82, re- 
spectively. We used eight sets of 2048 channel acousto-optical spectrometers (AOSs) that have a 
frequency resolution of 37kHz, corresponding to the velocity resolution of 0.48 km s~^. The typical 
system noise temperatures were in a range of 120-280 K during the observations. The achieved 
noise level in the spectrum was about 30 mK in T^, with a velocity resolution of 0.48 km s~^. The 
pointing accuracies were within 7" for the NH3 observations. 

The observational parameters are summarized in Tabled] for the 45 m telescope observations. 



2.2. Nobeyama Millimeter Array 



The aperture synthesis observations were carried out using the six-element NMA. We observed 
the H^^CO^(J = 1-0) line in the period from 2005 December to 2006 April. The phase track- 
ing center for our observations was set at the position of the IRAS source in B335. The field of 
view (i.e., the primary beam size) of the 10m dishes was about 79" at 87 GHz. Our observations 
were conducted with two array configurations (D and C configurations). For the maximum scale 
detectable in our observations, A/i?min ~ 71'.' l(i//87 GHz)~^(i?Tnin/1 0m)~^i the expected bright- 
ness recovered at the phase center is ~ 2.8% (jWilner &: Welchlll994l ) when a Gaussian brightness 
distribution having a FWHM of ^fwhm = ^/Bann is assumed. We used the SIS receivers as the 
front end and the digital spectro-correlator, New-FX, which has 1024 channels and bandwidth 
of 32 MHz (achieved velocity resolution was 0.108 km s~^), as the back end. We also obtained 
the contir iuum data simultaneo usly with the other spectro-correlator. Ultra- Wide Band Correlator 
(UWBC; lokumura et ahlboool ). which has 128 channels in the 1024MHz bandwidth mode. The 
system temperatures were typically in the range of 150-200 K in DSB. 

We used 3C345 and 3C454.3 as passband calibrators and B1923+210 as a gain calibrator. 
The absolute flux densities were calibrated using the flux density of B1923+210 determined by the 
bootstrapping method with Uranus and Neptune. The measured flux density of B1923+210 was 
1.26-1.61Jy at 87 GHz during the observation period. The uncertainty in our flux calibration was 



^ The Interactive Data Language 
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expected to be less than 10%. 

The NMA visibihty data were processed (cahbration, flagging, and continuum subtraction) 
using the UVPROCII package. For the continuum data, we merged the visibihties in each sideband 
obtained with the UWBC to construct an image. The cahbrated visibility data set was processed 
through imaging and image reconstruction with the MIRIAD package. The synthesized beam size 
is 5'.'5 X 4'.'3 (position angle (P-A.) = — 32?7) at 87 GHz with natural weighting. Imaging noise 
levels (1(t) of H^^CO^ and 87 GHz continuum data are 40mJybeam~^, with a velocity resolution 
of 0.108kms~^ and 0.56 mJy beam~^, respectively. 

The parameters for the NMA observations are summarized in Table [2j 



2.3. Combining the 45 m Telescope and the NMA Data 

We applied a data combining technique of interferometer and single-dish data in the Fourier 
domain (u-v domain) for our NMA and 45 m tel escope H^^CQ"|" data . We made the 45 m telescope 



image cube using a gridding kernel of Spheroidal (jSawada et al.ll2008l ) with a grid size of 7'.'92, which 



is the Nyquist spacing of a dish diameter of 45 m at 87 GHz. The velocity separation between 
the image channels is 0.108kms~^, which is the same as that of the NMA observations. The 
short spacing information in the u~v domain including the zero-baseline within the diameter of 
the primary antenna of NMA (= 10 m), which cannot inherently be obtained with the NMA 
observations, can be complemented with the 45 m telescope data. 



The fundamental theory and detailed algorithm of data combining are described in lKurono. Morita. &: Kamaza 



(|2009l ). In this method, we generated pseudo- visibilities from the Fourier transformed single-dish 
image data, which are deconvolved by a Gaussian beam with an FWHM of 25'.' 1 and multiplied by 
a primary beam of NMA approximated by a Gaussian function with an FWHM of 79". The single- 
dish visibilities that were generated were combined with the NMA data in the u-v domain. As a 
result, we can make synthesis images for a high spatial dynamic range and a high spatial resolution 
while recovering missing large-scale fluxes. We a.pplied the data optimizations to sen sitivities and 



relative weights between 45 m and NMA data (jKurono. Morita. &: Kamazaki I l2009l ) . We flnally 



obtained a combined image of B335 in the H^^CO"*" line emission with a synthesized beam size of 
5'.'6 X 4'.'4 (P.A. of 32?7) and with a total flux of 32.4 Jykms"^ that is comparable to that of the 
45 m telescope image, 30.0 Jykms"^. 

Furthermore, Y. Kurono et al. (2013, in preparation) describe the schemes that are applied to 
our observed data, and demonstrate the importance of relative (1) flux scaling, (2) sensitivity, and 
(3) weighting between interferometer and single-dish data in order to obtain reliable results. 
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3. RESULTS 

3.1. NH3 Line Emission 

The spectrum of the NH3(J, K) = (1, 1) transition, which consists of five hyperfine groups, and 
the main component of the (2, 2) transition were obtained with good signal-to-noise ratios for the 
45 m telescope. Figured] shows the spectra of the NH3(J, K) = (1, 1) and (2, 2) transitions obtained 
at the position of the IRAS source in B335. No emission in the (3,3) transition was detected. 

To estimate the optical depth and gas kinetic temperature, we analyzed the hyperfine struc- 
tures that are caused by the electric quadrupole moment of the nitrogen nucleus. We derived the 
peak main-beam temperatures (Tmb) and intrinsic velocity widths {Av) by fitting each line compo- 
nent with a Gaussian function. We assumed that all hyperfine components have equal beam filling 
factors and excitation temperatures. We estimated the optical depth of the NH3(1, 1) main compo- 
nent [T(l,l,m)], the rotational temperature [Trot (2, 2; 1, 1)] that describes the relative population 
between the (1, 1) and (2,2) levels, and the kinetic temperature (T^) by following the analysis 



in a previous work by iMangum et al.l ([1992) ■ With the parameters of the best-fit results (sum- 
marized in Table ED shown by the blue curves in Figure [TJ we obtained r(l,l,m) = 1.5 it 0.9, 
ri.ot(2,2; 1, 1) = 12.9 lb 1.1 K, and Tt = 15.1 it 1.5 K. The rotatio nal temperature is almost equal 



to the kinetic temperature for the range of Ti.ot(2, 2; 1, 1) < 20K (jDanbv et al.lll988l ). Therefore, 
we determined that the B335 core has a mean kinetic temperature of ~ 15 K over the beam area 
of the 45 m telescope (~ 80"), and used this temperature to estimate the isothermal sound speed 
and the column density from the molecular line data. 



3.2. 87GHz Continuum Emission 



We detected 87 GHz continuum emission with the NMA toward the IRAS source in B335, as 
shown in Figure E) The peak position was measured to be R.A. = 19''34™35!2, decl. = -h07°27^23.^8 
(B1950), which is con sistent with the peak position of the 1.3 mm continuum images (iHuard et al. 
I999I : lYen et al.ll2010l ). The 87 GHz continuum im age shows an elong ated emission structure per- 
pendicular to the outfiow lying from east to west (jHirano et al.lll988l l. Furthermore, a cavity- like 
distribution in the red-lobe side of the mole cular outfiow can be seen which is quite similar to 
that in the 1.3 mm continuum image shown bv lYen et al.l (|2010l ) and supports their suggestion that 
this distribution traces the wall of the outfiow cavity. The beam deconvolved size, peak intensity, 
and total flux density (above the 3a contour) are 9'.'2 x 6'.'4 (corresponding to 1380 AU x 960 AU), 
7.89 lb 0.56mJy beam~^, and 33.1 ± 0.17 mJy, respectively. 



As mentioned in IChandler &: Sargentl (|l993l ) , contribution from free-free emission is consid- 
ered to be negligible when taking into account the extrapolation from the flux at 3.6 cm. Thus, 
we determine that the 87 GHz continuum emission comes from a dust envelope surrounding a pro- 
tostellar object. Under the condition of being optically thin for thermal dust emission, the dust 
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envelope mass (Menv) was estimated using the equation Menv = Sy(fi / KyBy{T^), where Su is the 
total flux density, is the dust mass opacity coefficient, is the dust temperature, d is the 
source distance, and By is the Planck function. By combining our 87 GHz and image band 99 GH z 
measurements with th e flux densities at millime ter wavelengths estimated bv lKeene et alj (| 19831 ). 
Chandler et al.l ([l990|), and iHirano et al.l (|l992l ). we obtained a spectral index of a = 3.3 it 0.13 
for Sp oc v^. The spectral index gives /3 = 1.3 ± 0.13 for the emissivity law, k,j oc , using 
the a pproximated relation, /3 = a — 2, which i s valid for m illimeter wavelengths (IBeckwith et al. 
2000)- Thus, given K23ogHz = O.Olcm^g"^ (j Andrei 1 19941 ). we obtain a dust mass opacity of 
Kg? GHz = 2.8 X 10~^cm^g ~^. The dust envelope ma ss is estimated to be Menv ~ 0.19 Mq with a 
dust temperature of 30 K (jChandler &: Sargentlll993l ). 



3.3. H13QQ+ Lijjg Emission 

3.3.1. Integrated Intensity and Channel Maps 

The panels of Figure [3] show the integrated intensity images of the H^^CO+(J = 1-0) line 
emission over the LSR velocity range from 8.08 to 8.95 km s~^ made using the 45 m telescope data 
(left) and the combined 45 m telescope and NMA data (right). We should note that this combined 
map includes the effect of the primary beam attenuation of the NMA so that extended emission 
toward the outside of the image seen in the 45 m telescope image is not reproduced in the combined 
image. 

The H^^CO^ emission in the 45 m telescope map has a single-peaked spatial distribution and 
shows an elongation from north to south. The size of the H^^CO^ core above the 3(T level is 
~ 0.10 pc X 0.09 pc, with a P. A. of ~ 0°. The 45 m telescope plus NMA combined map shown in 
the right panel of Figure [3] clearly depicts the detailed structure of the inner core including large- 
scale flux distributions with a high resolution. The higher contours of > 12ct show an elongated 
distribution from north to south with a size of ~ 3000 AU x 1500 AU, which is believed to be an 
inner dense envelope associated with a central stellar source. The envelope has a double peak near 
the center, and the 87 GHz continuum source is located between the peaks. The elongation of the 



core and inner envelope is perpendicular to the molecular outflow axis (e.g.. iHirano et al.lll988l ). In 
the 45 m telescope map, there are faint ridges from the center along the P. A. of ~ 70°, ~ —45°, and 
— 135°, and they can be seen more clearly in the comb ined map. Taking the outflow direction 



and opening angle into account (~ 45°: IHirano et al.lll988l ). we believe that these ridges are related 
to the outflow activity. 

In the velocity channel maps obtained with the 45 m telescope shown in Figure IH the emission 
peaks are located at the north of the 87 GHz continuum source in the velocity range Vlsr = 8.19- 
8.62 km s~^, and at the south in the velocity range Vlsr = 8.73-8.84 km s~^. This velocity gradient 
is perpendicular to the outflow axis and can be interpreted as the rotational motion of the B335 
core. On the other hand, in the channel maps of combined data shown in Figure |5l it is difficult to 
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identify the corresponding velocity gradient because of the comphcated emission distribution. More 
detailed kinematics of the B335 core are discussed in Section 13.3.31 and 14.21 using position-velocity 
(PV) diagrams made from the 45 m telescope and combined images. 



3.3.2. Mass and Column Density Profile in the Core 

In order to estimate the column density of the B335 protostellar core, we analyzed the H^'^CO"'" 
images made from the 45-m telescope and the combined NMA plus 45 m telescope data. 

Under the local thermodynamic equilibrium (LTE) assumption, the H2 column density can be 
calculated using the following formula: 



JV„, 5.58 X 10'»^|jP3^ 



e-4.i6/Tex i_e-^ VKkms-i, 

where Tex is the excitation temperature, r is the optical depth of the H^'^CO^ line, and X(H^'^CO^) 
is the fractional abundance of H^'^CO' ^. When deriving of Equation ([1]), we used the perma- 



nent dipole moment = 4.07Debye (jHaese &: Wooda Il979l ) and the rotational constant B 



43377.17 MHz. To obtain a complementary expression for the combined synthesized image, we con- 
vert the antenna temperature into the flux density Sy in the Jy beam~^ through S^, = {2k-Q/X^^ Ob^lnb) 
where r^b is the beam solid angle given by 7r0maj^min/41n2. Thus, we have 

12 1 Tex 0.69 r 
iVH2 - 9.01 X 10^^ 



X(Hi3C0+) e-4.i6/Tex 1 _ e-^ 

9maj^min \ ^ / / Sydv \ ^^-2 

arcsec^ J \Jykms~^y 
We assumed the excitation temperature of 15 K from the NH3 analysis fSecti on 13.11) and the 



H^^C O+ line of the optically thin limit. We assumed X{Yi^^CO+) to be 8.3 x 10"^^ (IFrerking et al 



19871 ). The derived total mass from the 45 m telescope map (the left panel of Figured]) is ~ 1.2 Mq. 



We derived the radial column density profile in the B335 core from the column density map 
calculated by the above formula. The column densities were calculated from re-gridded images with 
cell sizes corresponding to the spatial resolutions of the images, i.e., 5" (750 AU) for the combined 
image and 25" (3750 AU) for the 45 m image. The column density profile was made as a function 
of the radius from the 87 GHz peak using the distributions of the estimated column densities over 
the re-gridded cells. In order to calculate the column density, the 45 m plus NMA image needed 
to be corrected for the primary beam attenuation, which increases the noise level in the outer 
region of the image. Hence, the column density profile in the outer region of the core was derived 
using the 45 m telescope data, and the 45 m plus NMA combined data were used to fill in the 
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inner region of < 3000 AU where the data cannot be obtained from 45m telescope data. Figure [6] 
represents the resulting H2 column density profiles obtained in the above procedure. As suggested 
in Section 13.3. H the B335 core could be affected by the outflow. To examine this effect we also 
made the column density profile by masking out the regions with P.A. of 67?5~112?5 and 247?5- 
292? 5. The maskin g angle of 45° was chosen to match the opening angle of the outflow: ^ 45° by 



Hirano et al.l (|l988l ) and ~ 41° by lHarvey et al.l (I2OOII ). Figure [6l^a) and (b) show the H2 column 



density profile without and with the masking, respectively. For both profiles, we can see that the 
column densities estimated from the 45 m and combined data are smoothly connected around the 
radius of ~ 3000 AU, and that the profiles in the inner radius are shallower than those in the outer 
radius. 

The column density profiles were fitted by two power-law functions of N{r) = NQ{r /rQ)~'^ , 
where s is the power-law index and A'^o is the column density at r = ro = 1500 AU. Since Figure [6] 
clearly shows two different slopes between the inner and outer regions, we estimated the turnover 
radius at which the power-law index changes. We evaluated the correlation coefficient of the power- 
law fitting in the range of rfit,i-rfit,0 5 where the inner radius (rgt^i) was variable and the outer radius 
('"fit,o) was fixed to rgt^o = 15,000 AU. As a result, we found that the correlation coefficient of the 
fitting decreased exceedingly when the fitting inner radius was taken inside of 4000 AU. Thus, we 
adopted a turnover radius of 4000 AU in this paper and performed the power-law fittings in the 
two regions of 750-4000 AU and 4000-15,000 AU separately 

For the inner region ranging from 750 to 4000 AU, we obtained N}i^{r) = (2.8 it 1.0) x 
1022(r/1500AU)-(°-^2±o.06) ^,^-2 without masking, and N^iir) = (2.6±1.0)xl022(r/1500 AU)-(°-^i=^°-''^) cm^^ 
with masking. For the outer region from 4000 to 15,000 AU, we obtained A'^h2('^) = (4-5 =b 1.3) x 
1022(r/1500AU)-(i-i°=^°-i^) cm-2 without masking and TVhz (r) = (4.0±1.3)xl022(r/1500 AU)~(i °i±°-^6) cm'^ 
with masking. The column density profile obtained with masking was estimated to be slightly shal- 
lower than that without the masks. From our fitting result, the column d ensity at r = 6000 AU is 



?a 9.7 X 10^^ cm-2, which is consistent with the estimate bv ISaito et all (|l999l ). A^HalSOOO AU) ~ 
6.3 X 10^^ cm~^. Using our results, we examined the density distribution of the B335 core (discussed 
in Section l4.1.ip and estimated the total mass of gas associated with the B335 core to be ^ 0.78Mq 
within the radius of 6000 AU. 



3.3.3. Position- Velocity Diagrams 

To examine the overall velocity structure in the B335 core, we made PV diagrams using the 
combined H^^CO"''(J = 1-0) line data as shown in Figure [71 with black contours and gray scales. 
We chose the two axes of P.A. = 0° (perpendicular to the outfiow axis) and P.A. = —90° (along the 
outfiow axis), passing through the 87 GHz continuum source. The comparison with the diagrams 
made from the low-spatial-resolution 45 m data overlaid with blue contours demonstrates a drastic 
effect of data combining. Along the P.A. of 0°, the 45 m telescope data show a slight velocity 
gradient, ~ 1.0kms~^pc~^ over a scale of 20,000 AU. This gradient is comparable to the typical 
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value of the cores in the Taurus dark cloud (jGoodman et al.lll993l ) . If the gradient represents a solid- 
body-like rotation of the core, then it is converted into an angular velocity of w ~ 3.2 x lO"^'^ rad s~^. 
On the other hand, the diagram of P.A. of —90° does not show a conspicuous velocity gradient. 

The diagrams of the 45 m and NMA combined data show line broadening and a double-peaked 
intensity profile at the core center (offset = 0"). With a broadening line width toward the core 
center, the diagram along the P.A. of 0° shows an overall symmetrical diamond-like structure (rep- 
resented by a dashed line in light gray) . We calculated intensity- weighted first and second moments 
defined by Mi = J I{v)vdv/ J I{v)dv and M2 = \J f I{v){v — MiYdv/ J I{v)dv, respectively. The 
second moment at the center (offset of 0") was M2 ~ 0.25 km s~^, whereas that at the offset of 20" 
was 0.17kms~^. The double-peaked profile at the core center is also slightly visible in the diagram 
made from the 45-m telescope data. This profile has asymmetrical features that are more enhanced 
in the blueshifted peak t han in the redshi f ted o ne. This blue-skewed profile is maybe due to the 
infall motion of the core ( Zhou et al.lll993l . 1 19941 ) because the H^^CO'''(J = 1 — 0) line emission is 
moderately optically thick around the systemic velocity (see Section I4.1.3P . 

The center of the overall symmetrical emission distribution (lowest contour), as represented 
with a dashed diamond in Figure [7|^ a) , is approximately ~ 8.4kms~^ and 0", which is considered 
to be the dynamical center of the B335 core. Since the LSR velocity of 8.4kms~^ is consistent 
with the peak velocities of the NH3 main components (Section 13. ip . we adopt an LSR velocity of 
8.4kms~^ as the systemic velocity of the protostellar system. 



4. DISCUSSION 

4.1. Density Structure of the Core 

4.1.1. Overall Structure of the Core 

When a molecular cloud core has a density distribution of p(r) oc r~P, the profile of the column 
density integrated along the line of sight has a relation of N{r) oc oc r^~P. Hence, our obtained 
indices of s ~ 0.5 in the inner part of the B335 core and of s ~ 1 in the outer part correspond to 
p ~ 1.5 and p ~ 2, respectively. We note that these indices barely varied even if we considered the 
influence of the cavity in the core made by the outflow. 

The column density profile derived from our 45 m and combined data with masking (Section 
I3.3.2P can be converted into the number density distribution: 

n(r) 1.2 X 10^ cm~^ 

(iooku)"'''' for 750 AU^r ^ 4000 AU, 
{mtm) for 4000 AU ^ r ^ 15,000 AU. 

Most previous studies have derived density profiles for the B335 star-forming core similar to 
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our observational results: inner and outer indices of —1.51 and —2.01, respectively, with a turnover 
radius of ~ 4000 AU. A power-law index of p ~ 1.95 for the density profile was estimated from 
single-dish observati ons of the C ^ ^0( J = 1-0) and H^^CO^(J = 1-0) line emissions (r adius range of 



4200- 25,000 AU) bv lSaito et al.l (|l999l ). Prom near-infrared extinction measurements. iHarvev et al 
(j200ll ) suggested that the B335 core has a constant power-law index of p = 1.91 over the region 
2000-15,000 AU, and otherwise displays insi de-out collapse with an infalling radius (see Section 



14. 2p of ~ 26" (corresponding to ~ 3900 AU). iHarvey et al.l (j2003al jbl) also suggested a single power 
law with p = 1.55 within the inner < 3000 AU region from interferometric observations of 1.2 
and 3.0 mm continuum. These results agree well with our estimate of the density profile from the 
H^ ^CO'^(J = 1- ) data . Meanwhile, detailed analysis of SCUBA 450 and 850 /im continuum maps 
bv IShirlev et all (j2002l ) showed that the model with a power-law index of p = 1.8 well describes 
the data as a best-fit result, and suggested a n infalling r a dius o f 1000 AU, which is smaller than 
estimates in our as well as previous studies. iDotv et al.l (j2010l ) recently conducted an unbiased 
fitting to the dust continuum observations toward B335. For the power-law density distribution, 
they obtained p = 1.5-1.9 throughout the envelope, although they did not find strong evidence of 
inside-out collapse with an infalling radius of ~ 1000 AU. 



4.1.2. Inner Structure of the Core 

We confirmed the above estimate of the density profile from our combined H^'^CO^(J = 1-0) 
line data using an analysis of the visibility data of the 87 GHz continuum emission observed with 
the NMA. This approach of directly examining in the u-v domain enables us to avoid possible 
artifacts caused by the deconvolution of interferometric images. 

When discussing the emission distribution of the inner region of an envelope observed with 
interferometers, we should consider the temperature distribution of the envelope. We consider 
the observed intensity as a function of the impact parameter, I,y(rb), from an optically thin dust 
envelope that has a spherical density distribution, p{r), and a dust temperature, 7d(r). If the 
density and temperature follow radial power laws, p oc r~'f and oc r~'^, and if we assume that 
the opacity does not vary along the line of sight, then the intensity profile also has a power-law 
profile, /i/(rb) oc r^^^^^'^ ^\ in the Rayleigh-Jeans regime. We assumed that the extent of source 
intensity /;/(rb) is sufficiently compact in the primary beam. The visibility amplitude as a function 
of u-v distance (6) can be given by V{b) oc 5^+^"^ for 3/2 < p + q < 3 (see the Appendix for 
details). 

Figure [8] shows a plot of binned visibility amplitudes as a function of the u-v distance for 
the 87 GHz continuum data. The binning is logarithmic, and the amplitudes are obtained as the 
vectorial average of the complex visibilities in each bin. The decrease in visibility amplitude with 
u-v distance can be interpreted as a power law. We fitted the visibility amplitude profile with a 
power-law function of the form A x b^^. The best-fit power-law index is = 1.13 it 0.06 (solid 
line in the plot), which can be converted into a power-law index with an intensity distribution of 
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p + q — 1 = 0.87. The power-law index of dust tempe rature distribution is expected to depend on 
the opacity, T^{r) oc L'^/^r"'?, where q = 2/(4 + /3) ( Dotv Sz LeuneJ Il994l ). Thus, the power-law 
index of the density distribution is estimated to be p = 1.49 for /3 = 1.3 (see Section [3.2p . which 
shows a good agreement with the estimate from the combined H^^CO"^ data (Section 14. l.ip . 



4.1.3. Uncertainties of Density Profile 

There are several uncertainties included in the derived density distribution. In this section, we 
examine the effect of opacity, assumption of uniform temperature, index conversion from column 
density to volume density profiles, and H^^CO^ fractional abundance. 

First, we investigate the effect of the optical depth of H^'^CO^ line emission, because the 
decrement of derived column de nsities in the inner radii might be due to larger optical depth within 



the 5" beam. ISaito et al.l (| 19991 ) indicated that the Hi3C0+ emission is optically thin even at the 
center of an 18" beam from the measured HC^^O"*" to H^'^CO^ peak intensity ratio of 5.8 and 
expected abundance ratio, X(H^^CO^)/X(HC^'^0^), of 5.5. Furthermore, the integrated intensity 
ratio of HC^'^O"*' to H^^C O^ implies that H^ ^CO"*" emission should be optically thin at the center. 



According to Figure 1 in ISaito et al.l (jl999l ). the line ratio actually ranges from ~ 2.8 to ~ 5.8 
across the velocity; the lowest line ratio indicates that H^^CO^ emission is marginally thick, r ~ 2 
at maximum in certain velocities. This opacity enhancement is limited only to a narrow velocity 
range, so it does not introduce a significant error when estimating the column density near the 
center, i.e., the underestimate. Therefore, the effect of the optical depth cannot account for the 
turnover around the radius of 4000 AU. 

Second, we assumed that the core has a uniform temperature distribution that does not have 
much effect on th e estimate of the column density. If we adopt the temperature profile from 



Evans et al.l (|2005l ). for example, then the difference in temperature between radii of 750 AU and 
20,000 AU should be a factor of ^ 2. From Equation ([2]), the uncertainty of the column density 
estimate expected by this temperature variation could be derived as it ~ 10% at most. 

Third, the derivation of the density profile of the core from the column density distribution 
depends on the radial finiteness of the core traced by the H^^CO"*" molecular line emission. For a 
spherical core with a size of i?out and a power-law density profile of p{r) = po{r/rQ)~P, the column 
density that is given by integrating the densities along the line of sight at the impact parameter 
can be expressed as 

/arccos(rb/-Rout) / n \ ^^^^ 

p{r)de = ATq - , (4) 

- arccos(rb/-Rout) \ ^0 / 

where 

/arccos(rb/-Rout) 
cosP9de, (5) 
- arccos(rb/-Rout) 

and the variable 6 is an angle between a radial vector r and the plane of the sky. Therefore, 
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the coefficient of A^o is also a function of r^, so that the column density profile deviates from a 
power-law of the form A^(rb) oc r^^^^ ■ Moreover, if the core has different power-law dependencies 
between the inner and outer regions, then the simple prediction of the column density profile in 
Equation @ is not exactly valid. Figure [9] shows the column density profiles of cores that have 
power-law density distributions with a cutoff radius -Rout • The dashed lines indicate the profiles for 
the predictions of oc r^'^^^. For the case of a single power-law density distribution of p{r) oc 
as shown in Figure a), the resulting column density coincides well with the dependency of oc r^^ 
at radii smaller than ~ 0.1i?out- However, it tends to depart from the power-law dependency and 
become steeper with increasing radius. Figure El^b) shows the column density profile for the core 
with density distributions of p{r) oc r~^'^ for r ^ 0.125i2out and p{r) oc for r > 0.125-Rout- This 
plot demonstrates that the column density profile of a core that has inner and outer regions with 
different power-law density distributions can be estimated to be steeper than the simple prediction 
of N{r\^) oc r^^'^^. We examined the uncertainties caused by these effects when converting the 
column density profile into the density profile, and we found that the true density profile is likely 
to be shallower by ~ 0.1 at most in the power-law index than for our estimates. 

Finally, as for the abundance of H^^CO"*", recent studies of chemical evolution in star-forming 
cores have shown that the abundance of HCO"*", which is a daughter spe c ies of CO, decrease s 



at radii that are smaller than the CO sublimation radius (jLee et al.l 12004 : lAikawa et al.l l2008l ). 



Nevertheless, it is shown that the abundance of HCO ^ at the inner radii i ncreases with the evolution 
of core collapse after a central stellar object is born. Evans et al.l (|2005l ) simulated a large number 
of molecular line profiles from B335 using various physical models. By the use of a self-consistent 
chemical model with core evolution, the result showed that the abundance of HCO"*" hardly changes 
along the core radius. The shallower profile at smaller radii can be affected by the radial distribution 
of the H^^CO"*" abundance, although it is not expected to be dominant. 

The possible effects of radial dependence in the H^'^CO^ abundance and the finite core radius 
oppose one another. It is difficult to estimate how they contribute to our data. Nevertheless, by 
taking into account the agreement of the column density profiles in the inner region of the B335 
core between the derivations from the combined H^'^CO"'' and dust continuum data, we determine 
that the effects of variation in optical depth and radius dependence in the H^^CO^ abundance 
are likely to be negligible. Moreover, the agreement between the profile indices derived from our 
data with those from the extinction analysis of near-infrared data indicates that the effect of 
these uncertainties, especially in the temperature and fractional abundance of H^'^CO"^, is not very 
considerable. In other words, our H^^CO"*" data well represent the column density structure of the 
core, which mostly ensures the validity of our analysis for the physics of star formation. 



4.2. Velocity Structure of the Core 



We performed model calculations of the PV diagrams and investigated which model reproduces 
the observed signatures of the PV diagram in the H^'^CO^ line well. For comparison with the 
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observed results, the model calculations were performed with different radial distributions of the 
infalling velocities and also for three different assumed central stellar masses. 



1.2.1. Model Calculations of PV Diagram 



We assume a contracting spherical star-forming core with rotation as described in detail be- 
low. Most of the parameters for the calculations simulating PV diagrams were taken from our 
observational results and previous studies. 

The core has a size of 20,000 AU and the radial density profile in Equation ^ shows power-law 
dependencies (p oc r~P) of p = 1.5 in the inner (r < 4000 AU) and p = 2 in the outer (r ^ 4000 AU) 
regions. Such a difference in power-law indices is naturally expected from the isothermal collapse 
model because of the boundary of the inner free-falling region and the outer region in which the 
condition at the stage of the formation of a central stellar object is expected to be conserved. The 
boundary radius is referred to as the infalling radius (nnf ), and therefore our derived density profile 
indicates rinf = 4000 AU. 

We also introduced the rotational motion of the core. The inclination angle of the rotation 
axis from the plane of the sky should be considered because it affects the line-of-sight velocity 
of rotational motion. We assumed that the rotati on axis correspond s to the outflow axis whos e 
inclination angle for B335 was estimated to be 10° (jHirano et al.lll988l ). 8 ±5° (ICabrit et al.lll988l ). 
and 9ibl° (iMoriarty-Schieven &: Snel]|ll989l ). We adopted an inclination angle of 10° for the rotation 
axis in our calculations. In the free-fall region with r ^ rinf, we adopted the increasing rotational 
velocity which is inversely proportional to the radius owing to the angular momentum conservation 
during collapse. Suc h a velocity field in infalling envelopes was indeed suggested observationally 
(e.g., L1551-IRS5 bv lMomose et al.lll998l ). The overall velocity gradient from the 45 m telescope 



observations, l.Okms" pc~ (w = 3.2 x 10~ rads~ ), can be regarded as the initial angular 
momentum of the core, which is still preserved in the outer radius of r ^ rjnf = 4000 AU. For the 
thermal line broadening, we set the line-of-sight velocity width of 0.15 km s^^ for H^^CO"^ molecules 
at Tk = 15 K. 

We calculated radiative transfer equations to estimate the relative intensity distributions and 
verify the effect of optical depth. Assuming a two-level state, the populations in rotational levels 
of an H^^CO"*" molecule were calculat e d usi ng the Einste i n A c oefficient and the collisional rate 
coefficients with H2 from iGerin et al.l (|2009l ) and iFloweii (|l999l ) . The final model PV diagrams 
were smoothed with Gaussian functions whose FWHMs are the same as actual resolutions of the 
combined image for both the spatial and velocity directions. 
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4-2.2. Comparison of Model Calculations of PV Diagrams with Observations 

We conducted simulations by varying the parameters of the inward velocity in the outer region 
{t ^ Hnf) of the core, f^^rinf, and the mass of the central stellar object to compare with observational 
results. We adopted the inward velocities in the outer region of fr^nnf = Okms~^, Cs, and 3.3 Cg. 
The centra l stell ar mass (Mc) for B335 was estimated in a wide range in previous studies, e.g., 
Choi et al.1 (Il99d ) estimated it to be 0.37 while O.O4M0 was suggested bv lVen et all (l2nid ). In 



our model calculations, we used three values of the central stellar mass: = O.OSMq, O.IMq, and 
O.15M0. Figure [To] shows the calculated PV diagrams perpendicular to the rotation axis passing 
through the core center. 

As for the inward velocities in the outer region of the core, the calculated PV diagrams for 
^f>nnf ~ Okms~^ well reproduce the features in the observed PV diagram. The inward velocity 
in the outer region affects the overall line width in the PV diagram. Although the PV diagrams 
calculated with Wr^rinf = Cs are also similar to the observed result, the overall velocity widths seem to 
be larger than that of the observed one. The simulated diagrams with Vr^n^i = 3.3 Cg undoubtedly 
disagree with the observation. We found that the central stellar mass mainly contributes line 
broadening at the core center, which means that the infalling motion is more dominant than the 
spin-up rotation owing to the conservation of angular momentum around the core center. For the 
diagrams with Vr^n^f — Okms~^, it seems that line wings at the center position are not enough in 
the case of Mc = O.OSMq while they are excessive in the case of Mc = O.ISMq. Consequently, we 
determine in our model calculations that the PV diagram with Vr^n^i — Okms~^ and Mc = O.IMq 
well represents the observed PV diagram. 

The calculated model PV diagram with Vr^^nf = Okms~^ and Mc = O.IMq shows that the 
peak is skewed toward the blueshifted velocity and has a shoulder in the profile around the systemic 
velocity at the center position. These features in the model diagram are due to the H^3co+ optical 
depth and are similar to the observational result. Hence, the suggestion discussed in Section [4.1.31 
that the H^^CO^(J = 1-0) line is marginally optically thick is probably reasonable. 

In addition, we can obtain the radial distribution of the specific angular momentum from the 
profile of rotation velocity in the model calculations of PV diagrams. Figure [IT] shows the radial 
distribution of specific angular momentum in the model that well represents the observation: the 
case of Vr^nnj = Okms"^ and Mc = O.IMq. We can see the spin-up rotation within the infall 
radius of rinf = 4000 AU as deviations from the solid-body rotation with an angular velocity of 
uj = 3.2 X 10~^^ rad s~^. In previous observational studies, the specific angular momenta of the core 
and e nvelope in B335 have be en measured at several radii from velocity gradients across the outflow 



axis. 



Yen et al.l ( 20ld . 1201 ll ) derived the specific angular momenta of < 7 x 10~^ kms~^ pc at a 
radius of 370 AU and 8 x 10~^ kms~^ pc at 90 AU from velocity gradients seen in the C^'^0(J = 2- 
1) and CS(J = 7-6) envelopes, respectiv ely. They discuss ed the evolution of specific angular 



momentum, including the measurements by lSaito et al.l ( 19991 ) of ~ 5.4 x 10 ^ km s ^ pc at 1000 AU 



and ~ 4.6 x 10 ^kms ^ pc at 20,000 AU. These measurements are also indicated in Figure [TT] 
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as open squares and open triangles. This plot shows that our model calculations successfully 
derived the radial distribution of the specific angular momentum consistent with measurements 
from previous observational studies. In particular, specific angular momenta that were recently 
estimated from high-resolution observations in submillimeter wavelengths at radii of 90 AU and 
370 AU show good agreement with the results of our model calculations. 



4.3. Comparison with Theoretical Models 



The density and velocity structures of a collapsing molecular cloud core are crucial to dis- 
tinguish theoretical models of gravitational collapse. We examine the observat ional results by 
comparing them with the key properties of the two sta r formation models, i.e., the IShul (jl977l ) and 



Larson-Penston solutions ( Larson 



196S 



Penstonlll969l . hereafter "LP solution" ) . Although the Shu 



and LP pictures of core contraction are similar in terms of the power-law dependencies of density 
distribution and inward motion, it is possible to distinguish between the two models by quantitative 
comparisons. 

The power-law dependency of the density profile in the B335 core qualitatively matches the 



similarity solutioii s for an isothermal spherical cloud in its post-protostar formatio n stage (IShu 



1977^ : lHunterill977l ). which has also been demonstrated by numerical simulations (e.g.. lOgino et al 



19991 ): i.e., p oc r~^/^ in the dynamical free-fall region and p oc in the outer region where 
the condition at the protostar formation stage should be conserved. The absolute density of the 
core is one of the key parameters used to discriminate between isothermal collapse models. We 
derived the number density of the B335 core as ~ 1.2 x 10^ cm~^ at the radius of 4000 AU. In 
Shu's inside-out picture, the density distribution in the outer region of the core is expected to 
correspond to the singular isothermal sphere, p{r) = c^/2iTGr~'^ . From this model, we obtain 
n(r = 4000 AU) = 9.2 x lO^cm"^ with T = 15 K. On the other hand, the "runaway collapse" of 
the LP solution during the core formation has a density that is 4.43 times higher than that in Shu's 
solution, so that we obtain 4.1 x 10^ cm~^ at r = 4000 AU. Our derived density is comparable to 



that predicted by Shu's solution with a factor of 
that of the LP solution. 



1.4, however, it is considerably smaller than 



The inward velocity in the outer region from our model calculations also supports Shu's solu- 
tion. As described above our model calculations of PV diagrams successfully explain the observed 
features within the uncertainties, in which the mass of the central stellar source and the inward 
velocity in the outer core region are estimated. Here, we focus on the inward motion in the outer 
region, which is one of the key characteristics to distinguish the gravitational collapse models. The 
static outer region of the core (fr^nnf = Okms~^) represents the velocity structure of Shu's solu- 
tion. On the other hand, numerical simulations of isothermal collapse models, which have been 
conducted to compare with similarity solutions, have showed that the inward velocity in the run- 
away collapse phase is not a constant but a function of age and radius. These inward velocities can 
be observed as intermediate pictures between those of the LP and Shu models, and are important 
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characteristics of the cohapse of an isothermal sphere. The infah velocity in t he runaway collaps e 
phase is determined by the initial ratio of gravitational force to pressure force (jOgino et al.lll999l ). 
The cloud core that begins to collapse from a marginally unstable isothermal gas sphere has a sub- 
sonic inward velocity in the outer region. Even from our model comparisons with the PV diagram 
of the observed data, we cannot exactly clarify if the outer region of the core has a subsonic inward 
velocity. However, it is plausible that there is an inward velocity > Cg in the outer region, such as 
Vr^nni = 3.3 Cs which represents LP solution, can be ruled out. 

Consequently, we suggest that the B335 core has initiated its gravitational collapse from a 
quasi-static initial condition similar to Shu's model. Otherwise, it is also possible to explain our 
observed results by means of a n isothermal collapse of a c loud core that has a mass slightly larger 



than the Bonner-Ebert mass (iFoster &: Chevalied 119931 : lOgino et al.lll999l ). In Shu's similarity 
solutions, an infalling radius provides a rough age of the cloud after a point source is formed at 
the core center, because the boundary between the infalling inner region and the outer region 
propagates outward as a rarefaction wave with a velocity of the isothermal sound speed. In our 
case, we take the turnover radius to be the infalling radius and obtain an age of ~ 8 x 10^ years 
which is comparable to the order of the age of the Class phase. Adopting the mass infall rate 
predicted by Shu's solution, M = 0.975Cs/G, along with the estimated age, we obtain a central 
stellar mass of ~ 0.2Mq. The estimate of the central stellar mass from our model calculation is 
closer to the prediction from Shu's solution than that from the LP model, which expects a 48 times 
higher mass infall rate. Therefore, the above comparisons between our results and the theoretical 
models indicate that the picture of Shu's model is more preferable than that of the LP solution for 
the B335 core. 



5. SUMMARY 

We presented a study of the dense molecular cloud core harboring the low-mass protostar, 
B335 in the H^^CO"'"(J =1 — 0) molecular line emission using the Nobeyama 45 m telescope and 
the NMA. Our main findings are summarized as follows. 

1. The single-dish observations revealed a dense core with a size of ~ 0.10 pc x 0.09 pc. Our 
analysis using a combining technique of single-dish and interferometer data revealed the structure 
of the inner dense envelope within the core with a high spatial resolution of ~ 750 AU. The envelope 
size is ~ 3000 AU x 1500 AU. Both of them have an elongated distribution toward the north-south 
direction, perpendicular to the outflow axis. The mass of the core is estimated to be ~ 1.2Mq 

2. We determined the radial column density profile of the B335 core and found a reliable 
difference between the power-law indices of the outer and inner regions of the dense core. The 
turnover radius is considered to be ~ 4000 AU, which is consistent with the infalling radius esti- 
mated in previous work. Our derived density profile, n(r) ~ 1.2 x 10^ (r/4000 AU)^^'^^ cm~^ for 
750 AU ^ r ^ 4000 AU and n{r) ^ 1.2 x 10^ (r/4000 AU)"^-°^ cm'^ for 4000 AU ^ r ^ 15,000 AU, 
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is better explained, both qualitatively and quantitatively, in the picture of Shu's self-similar solution 
than in that of the LP solution. 

3. The dense core shows a slight overall velocity gradient of ~ l.Okms"^ over the scale of 
20,000 AU across the outflow axis. This velocity gradient is considered to represent a solid-body 
rotation and corresponds to an angular velocity of a; ~ 3.2 x 10~^^rads~^. Our combined image 
also revealed detailed velocity structures in the dense core with a high resolution. The velocity 
structure of the B335 core can be well explained in terms of the collapse of an isothermal sphere, 
in which the core has an inner free-fall region and an outer region preserving the condition at the 
stage of protostar formation. 

4. We performed simple model calculations of PV diagrams to examine the observed diagrams. 
The model calculations successfully reproduce observational results, while suggesting a central 
stellar mass of ~ 0.1 Mq and a small inward velocity of ~ 0kms~^ in the outer region of the core 
> 4000 AU. 

5. Quantitative comparisons of density and velocity structures from the observational results 
with theoretical models show an agreement with Shu's quasi-static inside-out star formation. Fur- 
thermore, it is possible for the outer region of the B335 core to have a subsonic inward velocity. 
We concluded that a picture of Shu's solution or an isothermal collapse of a marginally stable 
Bonnor-Ebert sphere is suitable for the gravitational collapse of the B335 core. 

This study was based on observations at the Nobeyama Radio Observatory (NRO), which 
is a branch of the National Astronomical Observatory of Japan, National Institutes of Natural 
Sciences. The authors are grateful to the staff at the NRO for operating the NMA and the 45 
m telescope, helping us with the data reduction. We also thank an anonymous referee whose 
comments significantly improved the paper. Y.K. thanks Ken'ichi Tatematsu and the staff at the 
NAOJ Chile Observatory for their helpful comments, continuing interest, and encouragement. 

A. ANALYSIS FOR THE VISIBILITY FUNCTION OF DUST ENVELOPE 

In Section 14.1.21 we discuss the emission distribution of the inner region of the envelope using 
interferometric data in the u-v domain. The detailed expression of the analysis is described in this 
paper. 

For optically thin dust emission, the observed intensity from an envelope that has a spherical 
density distribution, p(r), and a dust temperature, 7d(r), as a function of the impact parameter, 
Tb, is written as 



where rout is the outer radius of the envelope. If the density and temperature follow radial de- 




dr, 



(Al) 
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pendencies of power-laws, p oc and oc r"'', and if we assume that the opacity does not 
vary along the line of sight, then the intensity also has a power-law profile, /(/(rt) oc r^^^'^'^ ^\ 
in the Rayleigh- Jeans regime. We assume that the intensity distribution, /,y(rb), is more compact 
in extent than the primary beam of interferometric observations, and gain variations during the 
observations are properly corrected. The visibility as a function of u-v distance, h = {y? + u^)^''^, 
can be given by the Hankel transform of the intensity distribution, 

V{h) = 2tt / 4(rb) Jo (27rrb6) r^dr^, (A2) 
Jo 

where Jq{z) is a zeroth-order Bessel function. 

Equation (|A2p is rewritten as a function of u-v distance 5, 

V{b)= Iu{rh) I exp[-27ri6rbCos(6i-a)]d(9rbdrb, (A3) 

where {l,m) = (rb cos ^, rb sin^) and [u^v] = (6 cos a, 6 sin a). By definition, a zeroth-order Bessel 
function is given by 

Jq{z) = — / exp (— cos 0) , (A4) 
2vr Jo 

so we obtain 

V{h) = 27r / /,(rb) Jo (2^rb6) r^dr^. (A5) 

JO 

This is the Hankel transform of the intensity distribution. We expect the intensity distribution to 
have a power-law dependency, this integral has a solution of the form (iGradshteyn h Ryzhiklll994l ) 

x^Jo(ax)dx = 2^a-^'-i— (A6) 

for ^ 

- 1< ;u < -, a > 0, (A7) 

where T{z) is the Gamma function: 



oo 



T{z)= / e-H'-'dt. (A8) 



Therefore, we obtain 

V{h) oc (A9) 

for 

-<p + g<3. (AlO) 

The visibilities of an intensity distribution with a spherically symmetric power law, / oc for 
1/2 < X < 2, are a power law in the u-v domain, V oc 
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Table 1. Summary of 45 m Telescope Observational Parameters 



Emission Line 


I/" 


Receiver 


fepBw'' 






^ e 


Mode* 


Area^ 




(GHz) 




(") 




(kms^^) 


(mK) 




(') 


NR3iJ,K) = (1,1)'' 


23.694506 


H22 


78 


0.82 


0.477 


31 


PS 


C 


NHsCJ.A-) = (2,2)'' 


23.722634 


H22 


78 


0.82 


0.476 


32 


PS 


C 


NH3(J,A:) = (3,3)'' 


23.870130 


H22 


78 


0.82 


0.473 


35 


PS 


c 


Hi3C0+(J = 1-0) 


86.75433 


BEARS 


18.5 


0.5 


0.108 


95 


OTP 


4x4 



''Rest frequency. 

''Half-power beam width for a Gaussian beam. 
"^Main-beam efficiency. 
''Velocity resolution. 

''Typical rms noise level of the spectrum. 

'Observing mode; PS denotes the position-switching observations and OTP denotes the On-The-Fly 
observing mode. 

^Size of the region for the mapping observations. "C" denotes the one-point observation toward the IRAS 
source at the core center. 

''Emission lines of three transitions were obtained simultaneously. 



Table 2. Summary of NMA Observational Parameters 



Emission Line i/* 


Configuration 


Phase Reference Center 


^pri'' 




Gain Calibrator 


Passband Calibrator 


(GHz) 


(u-^ Range (kA)) 


(B1950) 


(") 


(kms-1) 






h13cO+(J = 1-0) 86.75433 


D and C (2.89-47.2) 


19:34:35.1, 07:27:22.0 


78.9 


0.108 


B1923+210 


3C345, 3C454.3 



'^Rest frequency. 

''Primary beam size which is defined as full width at half-maximum for a circular Gaussian pattern. 
^Velocity resolution. 
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Table 3. Properties of NH3 Line Main Components 





NH3(1,1) " 






NH3(2,2) 




J- mb 




At; d 


mb 


WLSR " 


At; <i 


(K) 


(kms-1) 


(kms-i) 


(K) 


(kms-1) 


(kms-i) 


0.79 ± 0.06 


8.48 ± 0.05 


0.99 ± 0.08 


0.24 ± 0.07 


8.36 ± 0.08 


0.76 ±0.13 



'^Line properties for the brightest hyperfine components are shown. 
''Peak main-beam brightness temperature. 

■^LSR velocity at the peak brightness temperature by Gaussian fitting. 
^Velocity FWHM by Gaussian fitting. 




Fig. 1. — NH3(J, K) = (1,1) (upper panel) and {J,K) = (2,2) (lower panel) line profiles (black 
histogram) observed with the Nobeyama 45 m telescope toward the B335. Blue curves represent 
the best-fit results of Gaussian fitting (Section 13. ip . The rms noise levels are 42 mK in Tmb for both 
spectra. 
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RA offset (arcsec) 

Fig. 2. — 87 GHz continuum image toward the IRAS source in B335 observed with the NMA. 
Contours start from the 3a level with a 3a interval, where la = 0.56mJy beam~^. Dashed contours 
show the —3a level. Solid and dashed arrows in the image indicate the outflow directions of blue 
and red lobes, respectively. The ellipse at the bottom right corner indicates the synthesized beam 
size (5'.'33 x 4'.'31 with a position angle of -37?2). 
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Fig. 3.— Total integrated intensity maps of B335 in the H^^CO+(J = 1-0) line emission obtained 
with the 45 m telescope (left) and by combining the 45 m telescope and NMA data (right). Contour 
intervals are 3o", starting from the 3a levels of the images where la noise levels are 62mKkms~^ 
in Tjiib for the 45 m telescope map and 24mJy beam^^ kms^^ for the combined map. Dashed 
contours show the —3a level. The open star in each map is the peak position of the 87 GHz 
continuum emission observed with the NMA. The beam size for each map is shown as a filled circle 
or filled ellipse at the bottom right corner. Dotted circle in the combined image indicates the field 
of view, i.e., FWHM primary beam size of the NMA observations. 
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Fig. 4. — Velocity channel maps of B335 in the H^'^CO^(J = 1-0) line emission taken with the 45 
m telescope. The central LSR velocity in kms~^ for each channel is shown in the upper left corner. 
Contour intervals are the same as in Figure El and the la noise level is 0.19 K in Tmb- The open 
star, filled circle in each panel, and dotted circle are the same as in Figure El 
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Fig. 5. — Velocity channel maps (CLEANed synthesized images) of B335 in the H^'^CO^(J = 1-0) 
line emission obtained by combining the 45 m telescope and NMA data. The central LSR velocity 
in kms~^ for each channel is shown in the upper left corner. Contour intervals are the same as 
in Figure [3] and la noise level is 45mJybeam^^. The open star, filled ellipse in each panel, and 
dotted circle are the same as in Figure [3l 
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Fig. 6. — Radial profile of H2 column density in the B335 core made from the H^'^CO^(J = 1-0) 
image obtained with the 45 m telescope only (filled circles) and the combined 45 ni telescope and 
NMA data (filled triangles). Panel (a) shows the profile without the mask in the cavity regions 
and panel (b) shows that with the mask. Dot-dashed and dashed lines in the plot indicate the 
best-fit results of power-law fitting with the fitting range of 750-4000 AU and 4000-15,000 AU, 
respectively. Dotted curves indicate approximated beam patterns of 45 m telescope image and the 
combined image, Gaussian functions with the FWHM of 25'.'0 and 5''0, respectively. 



- 33 - 




Fig. 7. — Position- velocity diagrams of the H^^CO^(J = 1-0) line emission for the B335 core made 
from the 45 m telescope image data (blue contours) and the combined data (gray scale with black 
contours). Panel (a) and (b), respectively, show the diagrams along P. A. = 0° (perpendicular to 
the outflow axis) and P. A. = 90° (along the outflow axis), passing through the 87 GHz continuum 
source. 
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Fig. 8. — Binned visibility amplitude of the 87 GHz continuum data taken with the NMA as a 
function of u-v distance. Binning is logarithmic and amplitudes are obtained as vectorial average 
of the complex visibilities in each bin. A solid line shows the best-fit power-law function. 
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Fig. 9. — Column density profiles of core having power-law density distributions with a cutoff 

radius of .Rout- Panel (a) shows the case of a density distribution of p(r) oc and panel (b) shows 
the case of a density distribution of p{r) oc r^^-^ (r ^ 0.125i?out) and p(r) oc (r > 0.125-Rout)- 
Dashed lines indicate column density profiles expected from simple predictions of oc for (a), 
and oc r'^-^ {r ^ 0.125i?out) and oc r'^ (r > 0.125i2out) for (b). 
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Fig. 10. — Panels (a)-(i) show simulated PV diagrams by the simple calculations using models 
of contracting core with rotational motion. The cutting axis is perpendicular to the rotation axis. 
Panels (a), (b), and (c) show the model PV diagrams with Vr^^^f = Okms"-^, Cg, and 3.3 Cg with 
the central stellar mass of Mc = O.OSM©, respectively. Panels (d)-(f) show the diagrams in the 
cases of Mc = O.IMq and panels (g)-(i) show the cases of Mc = O.ISM©. PV diagram of the 
H^^CO"'"( J = 1-0) line emission from B335 perpendicular to the outflow axis is shown in panel (j) 
with black contours. Contour intervals of model PV diagrams relative to the peaks are the same 
as those for the observed diagram. 
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Fig. 11. — Radial distribution of specific angular momentum in our model calculation (Section 
14. 2p for the case of Vr-^n^i — Okms~^ and Mc = O.IMq (cross marks with solid line) as a function 
of radius. Open squares represent specific angular mor nenta measured f r om ve locity gradients seen 
in the C^^O{J = 2-1) and CS(J = 7-6) envelopes bv Iven et aP (l2ninl . boilh . and open triangles 
denote spec ific angular mome nta measured from C^^O(J = 1-0) core and H^^CO"'"(J = 1-0) 
envelope bv lSaito et al.l (|l999l ). A dashed line indicates the profile of specific angular momentum 
of a; = 3.2 x lO-i'^rads'^ 



